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Abstract 

We show that correlation functions for branched polymers correspond to those for 3 
theory with a single mass insertion, not those for the (ft 3 theory themselves, as has been 
widely believed. In particular, the two-point function behaves as 1/p 4 , not as l/p 2 - 
This behavior is consistent with the fact that the Hausdorff dimension of the branched 
polymer is four. 
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1 Introduction 



Branched polymers are the simplest generalization of the random walk and have been studied 
extensively. It is of great importance not only in statistical physics but also in particle 
physics, in particular, for understanding the critical behavior of random surface models and 
quantum gravity (see, for example, U). In the paper 0, we have studied the dynamics 
of the type IIB matrix model of a constructive formulation of superstring (0 and see [|J 
for review). In our matrix model approach, the eigenvalues of matrices are interpreted as 
space-time coordinates. In these investigations, we find the system of branched polymers in 
a simple approximation. Although it is far from the flat four dimensional manifold, branched 
polymers share the same (fractal) dimensionality four with our space-time. It might be the 
first indication that superstring can explain the dimensionality of our space-time. 

In this letter, we comment on a field theoretic description of branched polymers. It is 
well-known that a system of random walks is described by a free scalar field theory if there 
is no effect of self-avoidance. Similarly it is widely believed that the system of branched 
polymers becomes a scalar field theory with a three-point coupling, that is, 3 scalar field 
theory (see, for example, [l|]). We will, however, show that it is not so by treating the 
universal part of the partition function carefully. The system of branched polymers without 
self-avoidance can be exactly solvable by introducing the grand canonical ensemble and using 
the so called Schwinger-Dyson technique. In order to extract the correct large N limit (N 
is the system size) or the thermodynamic limit, we have to check that the grand canonical 
ensemble is dominated by the larger size system. In other words, we have to extract the 
universal part. 

On the contrary to a naive summation which makes us conclude that the system of 
branched polymers is described by a 3 scalar field theory, a careful treatment shows that 
we need a single mass insertion in each n-point function of the 3 scalar field theory. Mass 
insertion here means a change of a propagator in each n-point function from an ordinary one, 
1/p 2 , to 1/p 4 . In particular, the two-point function behaves as 1/p 4 , not l/p 2 . Let us count 
the number of points which lie within distance R from a fixed point in D(> 4) dimensions. 
In the random walk, it can be estimated as R 2 by using the two point function. We obtain 
R A for branched polymers by using the 1/p 4 type propagator. So our finding are consistent 
with the claim that branched polymers are four dimensional fractals. A multi-point Green 
function is given by a sum of graphs of the corresponding Green function for <ft 3 scalar field 
theory at tree level with a single mass insertion in each graph. 

The organization of this paper is the following. First we define a canonical ensemble 
for a system of branched polymers (sec. 2.1) and then introduce grand canonical ensembles 
(sec. 2.2). We emphasize that a definition of grand canonical ensemble is not unique. In 
sec. 2.3, we solve Schwinger-Dyson equations for the grand canonical ensembles and obtain 
naive results of the correlation functions of a scalar 3 theory. In sec. 2.4, we consider the 
thermodynamic limit and give the correct universal behavior of the correlation functions. 
Finally in section 3, we conclude and give an interpretation of our result. 

2 Branched polymer dynamics 
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2.1 Canonical Ensemble 



Branched polymers are a statistical system of N points connected by N — 1 bonds whose 
lengths are of order clq. The canonical partition function is defined as 



I r N ■ 

zn = W] e /n^ n /(*• 

G:tree graph i=l (ij):bonf of G 



where f(x) is a function assigned to each bond in each graph, and it damps sufficiently fast 
at long distance compared to the typical length scale ao- (See, for example, figure 1.) The 
presence of the factor 1/N\ is due to the fact that the N points are regarded identical. 

We can calculate a partition function for each N, by counting the number of all possible 
tree graphs: 
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where f(p) is a Fourier transform of f{x) 

f{p) = I dx e*~ f(x), 



(6) 
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and V is the total volume of the system. A derivation of the general form (|6|) is given in the 
appendix. 

We define an (unnormalized) m-point correlation function of density operators as 
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Figure 1: Form of the function f(x). It damps sufficiently fast at long distance compared to 
the typical length scale a . 
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graph i=l (ij) 

where the density operator is defined by 

N 



p(x) = ^5 (d) (^-^)- (9) 



i=l 



2.2 Grand canonical ensemble 

We then define partition functions and m-point correlation functions in the grand canonical 
ensemble in the following equations: 

oo 

z^i = Y,n 1 $z n , (10) 

N=l 

oo 

GZ,(x\---,x m ) = ^N 1 k»G%(x\---,x™), (11) 

N=l 

where k is fugacity. There exist various ways of defining grand canonical ensemble cor- 
responding to the freedom to choose different weights for each fixed N sector. Here, we 
assigned an extra iV-dependent factor, N l , in addition to the usual one, Kq. The criterion 
for a "good" grand canonical ensemble is such that we can take the correct thermodynamic 
limit, or, in other words, we can correctly take the universal part in the sum (|T0|) or (|TT|) . 
That is, the correlation functions in the grand canonical ensemble at the critical value of 
fugacity should reproduce those in the canonical ensemble for large N: 

r<m (~~<m 

lim )fLH = lim -J*l. (12) 

Af^oo Zn k o— «o, c Z Kq i 

We illustrate the above mentioned criterion by taking partition function as an example. 
Since the canonical ensemble partition function (|6]) behaves at large N as 

"n ~ -^^fiOf-'V, (13) 



the grand canonical partition function is approximated by 



V2nf(0) Jo 
where 



V2irj[0) n=1 

oo 

dN N l ~ 5 ' 2 e ~ NAK/K % (15) 



V 



k = /(0)«o, (16) 
« c = e~\ (17) 
An = n r — K. (18) 
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If we take I sufficiently large, the integrand in eq. ([D]) has a peak at iV ~ k c /Ak and we 
can make N large by approaching k to k c . On the other hand, if / is not sufficiently large, a 
non-universal small N behavior dominates the summation and we cannot obtain the correct 
answer of the large N limit by a grand canonical ensemble. 



2.3 Schwinger-Dyson eq. 

In this subsection, we recapitulate the (naive) arguments that the correlation functions 
for branched polymers are given by massless <f) 3 theory JT|[@. Let's consider the following 
correlation functions, which are suitable for a Schwinger-Dyson analysis: 

oo N . N 

F m (x\---,x m ) = Era E JUdy* 5to(y l -J)---5W{y m -x m ) 

N=l V /" G:trce graph i=l 

x n f(y l -y 3 )- (19) 

(ij'):bond of G 

The factor ^ N ^ m y means that (N — m) points other than the fixed m points are regarded 
identical. We can see from eq. flTTD and eq. (§) that 

F m (x\---,x m )~G™ 0jl=0 (x\--- } x m ) (20) 

in the large N limit. We write a Fourier transform of F m (x 1 , ■ ■ • , x m ) as F m (p l , ■ ■ ■ ,p m_1 ): 

(27r) d 6 w (p 1 + ■■■+ p m )F m {p\ ■ ■ ■ ,p m ~ 1 ) 

= J d d x l ■ ■ ■ d d x m e iplxl ■ ■ ■ e ipmxm F m (x 1 , ■ ■ ■ , x m ). (21) 

Schwinger Dyson equation for 1-point function F l becomes 

b = ne b , (22) 

where 

b = f(p)F\ (23) 

as can be seen from figure ||. Figure | illustrates eq. ( p2|) . 
At the critical point, 

b c = 1, (24) 
k c = e-\ (25) 

1-point function db/dk diverges. Near this critical point, iV becomes large; 

Ab ~ \Plt\[K^ ~ 1 / VN, (26) 

where 

Ab = b c - b, (27) 
Ak = k c - k. (28) 
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Figure 2: Schwinger- Dyson equation for one point function. A grey blob and a black point 
mean Fi and kq, respectively. 




Figure 3: Schwinger-Dyson equation, k = be b . At the critical point, b c = 1, k c = e l . 



Next, we consider the 2-point function F 2 (p). When we pick up any two points on a tree 
graph, we can fix the path connecting these two points. Thus, as can be seen from figure [|, 
2-point function is calculated to be 

oo 

F 2 (p) = E/(P) S (^) S+1 (29) 



8=1 



b 2 h(p) 



f(0)(l-bh(p)) 



(30) 



where 



h(p) = /(p)//(0) (31) 

= l-ca 2 p 2 + ■■■. (32) 

Here c is a positive constant of order one. Recall that f(x) damps rapidly out of the 
region < x < a , as shown in figure 1. 

Near the critical point, b ~ b c = 1, 2-point correlation function behaves as 

F 2 {p) ~ \ i ~ 90 (33) 

v ' Ab + ca^p 2 v ; 
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Figure 4: 2-point function F 2 is made out of 1-point function F 1 , which is written by a gray 
blob. 



1 

AT- 1/2 + m 2 p 2- ( 34 ) 

Here we used eq. (|26|) . Thus, the correlation length is £ = clqN 1 I a , which implies the 
Hausdorff dimension of branched polymer is four. If the relevant length scale is shorter than 
the correlation length, F 2 {p) ~ 1/p 2 behaves like a propagator of a massless scalar field. Let 
us consider the following region: 

a < x < f = aoN 1 ' 4 . (35) 

cio gives an ultraviolet cut-off whereas £ gives an infrared cut-off length over which correlation 
functions damp rapidly. Note that, in this region and near the critical point, the following 
inequality holds: 

1 « i^sm « rb- < 36 > 

Finally, we consider correlation functions of m > 2. As in the case of the two-point 
function, m points are fixed on each tree graph. We can uniquely fix the path connecting 
them on each graph. Therefore, an m-point function F m is represented as a summation over 
all tree diagrams with m fixed points in which F 2 appear as propagators. For example, 

F 3 (p,q) = {b/f{0))- 2 F 2 { P )F 2 { q )F 2 { P + q) 

+(b/f(0))- 1 [F 2 {p)F 2 (q) + F 2 (p)F 2 (p + q) + F 2 (q)F 2 (p + q)}, (37) 



as we can see from figure g. However, because of the inequality (|3q), the diagrams with the 
maximum number of propagators dominate. In the case of 3-point function, the first term 




+ 



Figure 5: 3-point function F 3 is a summation of tree diagrams whose propagators are F 2 
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in eq. ( 57 ) dominates. In general, we obtain a naive result for m-point correlation functions; 

F m ~ correlation functions of massless 3 theory at tree level. (38) 

In the next subsection, we will see that these naive results do not correspond to the correct 
thermodynamic results 

2.4 Correlation functions in thermodynamic limit 

As we mentioned in eq. (|20|), G™ l=0 is equal to F m in the large N limit: 

<W=o(p\ • • • >P m_1 ) - ^(P^ ■ ■ ■^ m ~ 1 )- (39) 

Then, let us consider the m-point correlation functions with I > 1. ^From the definition 
(PI), they can be obtained by applying l-th derivative to the I = case: 

GZM'-'-'P™' 1 ) = {^) l G2,i= Q (jp\---,P m - X ) (40) 

= ^M^Up 1 ^---^ 1 ) ( 41 ) 

2-point function with / = 1 is given by 



1 



(1 - b)(l - bh(p)) 2 
1 1 



(1-6) p 4 



(43) 
(44) 



Here, we have used the inequality (^) since we are interested in the correlation functions 
in the physical region (|3^) near the critical point. This behavior is different from that of 
G 2 kq1=q ~ l/p 2 - Similarly for I > 1, the behavior of 2-point function becomes 

d '-^ {p) ~ (1 - 6)*-' (1 - bh(p)f (45) 

1 1 (46) 



(1 - b) 21 - 1 J9 4 

and coincides with the I = 1 result. Due to the inequality (|36|), the derivative is dominated 
to act on not on x _ hh ^ ■ Thus, their p-dependences are all the same. This is the correct 
thermodynamic limit. That is, we should consider a grand canonical correlation function 
with I > 1, otherwise a non- universal small N behavior affects the summation and we cannot 
obtain the universal result. Note that the behavior of G 2 Kol>1 (p) ~ l/p 4 is consistent with 
the fact that the Hausdorff dimension of the branched polymer is four. A naive argument 
expected from the figure ^ is that the effect of branching can be absorbed by renormalizing 
the mass. If so, the propagator behaves as that of random walks. We discuss in section 3 
why this argument is not correct. 
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Similarly, 3-point functions become 

G 3 K0tl=0 ( P ,q) ~ g(p)g(q)g(p + q), (47) 

G 3 K0 ,i>i(p, q) ~ g(p)'g(q)g(p + q)+ g(p)g(q)'g(p + q) + g(p)g(q)g(p + q)\ (48) 

where 

9(p) = rdsw~?" (49) 

9(p) ' = (i - ikcp)) 2 ~ (50) 

Note that only one propagator in a graph is replaced by g'(p), since the derivative ^ is 
dominated to act on the factor y^, as in the case of 2-point functions. Therefore, for 
m-point correlation functions, 

G™ l=Q ~ correlation functions for 3 theory at tree level (51) 
G™ 1>\ ~ correlation functions for 3 theory at tree level with a mass insertion (52) 



and the correct correlation function in the thermodynamic limit should be given by eq. (|52|) , 
not by eq. fl5jj). 

As a consistency check, the following relation between an (m+l)-point function and an 
m-point function must hold: 

G m+1 (p 1 ,---,p m - 1 ,p m = 0) =G m {p 1 ,---,p m - 1 ). (53) 

It actually holds because in the L.H.S. of eq. (0), the special class of diagrams dominate in 
which the m-th end point is attached to the propagator g'{p m ). It is due to the inequality 
g(p = 0) g(p = 0)'. Then, it is equal to the R.H.S. of eq. (|53"D . 



3 Conclusion and discussion 

In this letter we have shown that the correlation functions for branched polymers are given 
by those for 3 theory at tree level with a single mass insertion if we correctly take the 
thermodynamic limit. It is not given by those for 3 theory at tree level themselves, as has 
been widely believed. 

Our result can be interpreted as follows. Since the 1-point function behaves as 

Gjv ~ N-W(ef(0)) N (54) 

at large N (see eq. (|T5|)), we obtain a relation; 

G\j 3> J dN 1 G l N i G^ N _ N /y (55) 

This relation reminds us of the situation in the two-dimensional quantum gravity (see, for 
example, H). Two dimensional quantum gravity is known to describe a fractal space-time. 
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It consists of numerous tiny baby universes and a single mother universe. In our case, this 
analogy means that if we divide any tree graph into two by cutting a bond, we find only finite 
points in one of them and the most of the points belong to the other. In the case of m-point 
functions, there is a single mother universe on a blob in the path connecting these m points. 
Let's consider the case of the two-point function as a simplest example. Naively the figure [| 
implies that the propagator behaves as that of the simpler random walk if each blob is on an 
equal footing. However one of the blobs in figure f| becomes the mother universe consisting 
of infinitely many points and it is entirely different from the other blobs. Therefore we have 
to divide the path into two parts by cutting out the mother universe. Each part can be 
considered as a propagator of random walks whose weights are dressed by blobs. This is 
the reason why we have obtained the propagator behaving as 1/p 4 instead of 1/p 2 . In other 
words, the mother universe (the blob with infinitely many points) corresponds to a mass 
insertion and the other blobs renormalize the correlation length £ of random walks from 
a iV 1/2 to a N l /\ We can also apply a similar argument for the higher-point functions. One 
of the blobs becomes the mother universe, which corresponds to the single mass insertion in 
the (ft 3 scalar field theory. 

Acknowledgments 

We would like to thank Dr. S. Higuchi for correspondence. This work is supported in part 
by the Grant-in-Aid for Scientific Research from the Ministry of Education, Science and 
Culture of Japan and by the National Science Foundation under Grant No. PHYS94-07194. 



Appendix 



In this appendix, we derive the canonical-ensemble partition function from the Schwinger- 
Dyson equation (^). Let us solve b as a form of expansion in k. Each coefficient is calculated 
to be 

— / dn-^— = — I db(l - b)e~ b - — (56) 



iV 



N-l 



Hence, 



From the definitions of ([!]), fll9D and (p3|) , b is expanded as 



Comparing these two expansions, we get the result of eq. 



(57) 



oo AfN— 1 

» = E V K " < 58) 

JV=1 iV ' 



00 N 

b = /(0) E 77 Z » K o- (59) 
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